Signatures of glass formation in a fluidized bed of hard spheres 
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We demonstrate that a fluidized bed of hard spheres during defluidization displays properties 
associated with formation of a glass. The final state is rate dependent, and as this state is ap- 
proached, the bed exhibits heterogeneity with increasing time and length scales. The formation of 
a glass results in the loss of fluidization and an arrest of macroscopic particle motion. Microscopic 
motion persists in this state, but the bed can be jammed by application of a small increase in flow 
rate. Thus a fluidized bed can serve as a test system for studies of glass formation and jamming. 

PACS numbers: 64.70.Pf,81.05.Rm,47.55.Kf,45.70.-n 



Despite a century of study, there is no unified theory 
for the formation of a glass state from a liquid cooled 
below its freezing point (supercooled) Q, However, 
some signatures of the glass transition have been identi- 
fied. One is the dependence of the volume of the glass on 
the cooling rate: the slower a liquid is cooled, the greater 
density it achieves when it forms a glass. Another sig- 
nature is a rapid increase in a characteristic relaxation 
time (or viscosity) of the liquid during cooling as the 
glass transition temperature T g is approached. In col- 
loidal and hard sphere systems the transition is reached 
by increasing volume fraction 4> rather than decreasing 
temperature; the transition occurs at (j) g ss 0.58 0,0,0- 

In liquids El H, B, E|, colloids 0, E3, and simula- 
tions it has been observed that as the glass 
transition is approached, the system dynamics becomes 
increasingly heterogeneous. Both the size of spatially cor- 
related regions and the time scale for rearrangement of 
these regions increase rapidly near the glass transition. 

We examine a fluidized bed 0, 0] , which is a verti- 
cal column of particles with an upward flow of fluid that 
maintains the particles in motion. As the fluid flow rate 
Q is decreased, this non-equilibrium system exhibits 
characteristic signatures of the equilibrium glass transi- 
tion. We find that two different volume fractions are 
important in the loss of fluidization: (f> g = 0.582, where 
bulk quantities change, and 4> a = 0.594, where <f> becomes 
nearly independent of Q and only small scale motion per- 
sists; this motion can be stopped by a slight increase in 
flow rate, jamming the system. Our observed (f> g is close 
to the volume fraction at which a glass transition is ob- 
served in other hard sphere systems 0. Q|. 0] . The 
observed 4>a is close to the value obtained for slow deflu- 
idization of particles in a gas fluidized bed 0] , and close 
to the volume fraction where the system has the maximal 
number of statistically independent regions . 

Experiment — Water flows upward at a volume flow 
rate Q through a vertical column of glass spheres in 
a square bore glass tube of cross sectional area A — 
5.81 cm 2 . There are 4 x 10 6 glass spheres with diameter 
d = 250 ± 8 fim and density p p — 2.47 g/cm 3 . Flow rate 
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FIG. 1: The volume fraction achieved during defluidization 
depends on the ramp rate (cf. values next to the curves), just 
as in the glass transition in colloids and in hard sphere models. 
In the top two curves the ramping is so rapid that the final 
state is achieved through sedimentation. Inset: magnification 
of the region near the arrest transition at <f) a = 0.594; only 
the three slowest ramp rates are shown. 



fluctuations are smaller than 0.3%. To obtain uniform 
flow, fluid passes into the bottom of the column through 
a nylon mesh (5 /J,m weave, open area 0.75%, Nitex mesh, 
Sefar America) . The average height h of the bed top sur- 
face above the distributor is measured to determine the 
average bed solid volume fraction, <f> = M p / p p Ah, where 
M p is the total mass of the particles. We measure the 
fluid pressure drop AP from the bottom to the top of 
the bed; values are normalized by the buoyant weight 
of the grains. The dynamics of the bed are studied by 
imaging the side of the bed and by using Diffusing Wave 
Spectroscopy (DWS) [2(J to probe the interior of the bed. 

Signatures of glass formation — The dependence of <p 
on the ramping rate (Fig. ^) is similar to the dependence 
of glass volume on cooling rate in supercooled liquids that 
form glasses In a fluidized bed for sufficiently large 
Q, the grains are mobile (2^ and AP = 1. However, as 
Q decreases, particle mobility decreases and <p increases 
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FIG. 2: (a) Dynamical heterogeneity: as the glass state is 
approached, transiently immobile regions (black) and tran- 
siently mobile regions (white) grow in size. Each panel is the 
difference between two images (1.5 x 1.5 cm 2 ) of the side of 
the bed taken AT = 0.25 s apart, (b) Raster scans for the 
row of pixels indicated by an arrow in (a) show that the life- 
times of the correlated regions increases as the transition is 
approached. The flow rates for (f> — 0.568 and cf> = 0.590 were 
respectively 33.02 and 26.19 mL/min. In this figure and in 
Fig. 3 the data were obtained at (f> values reached by decreas- 
ing Q at the slowest ramp rate in Fig. 



until the defluidization transition is reached at cj) a (the 
knee of the curve in the inset of Fig. 1); (f> a decreases 
only 1% as the "cooling rate" dQ/dt is decreased by more 
than two orders of magnitude. 

Another signature of glass formation is dynamical het- 
erogeneity. As in other glass-forming systems 0, 0, El 

mlli, the 

size of the mobile and immobile regions in 
the fluidized bed increases as increases 0, as shown 
by the images in Fig. Eta) and by the correlation range 
£ in Fig. EH a), where £ is given by the 1/e point of an 
azimuthal average of 2D autocorrelation of difference im- 
ages separated by time AT. The heterogeneities observed 
both at the top and the side of the bed are similar to vol- 
canos and channels observed in gas fluidized beds |23| . 

The decrease in £ near </> 9 is similar to the decrease 
in length scale of cooperative regions observed in colloid 
experiments 01 an d lattice model simulations 01 ■ The 
decrease in £ that we observe arises because the dynamics 
has slowed enough near <j) g so that no discernible motion 
occurs during the measurement time AT. Motion can 
be observed for increased AT, but th en £ remains close 
to zero until (j) « (j) g (cf. Fig. |3Ja) [21j); most parti- 
cles in the bed have moved during this AT and thus the 
correlation measurement picks up only the minimum ob- 
servable motion. The lifetime r of mobile and immobile 
regions increases as 4> increases, as shown in Fig. Efb) 
and in Fig. |3Jb), where r is obtained by averaging over 
all pixels (i,j) (for each pixel the time Tij is given by the 
1/e point of the time autocorrelation function of the in- 
tensity). The standard deviation of the distribution of 
decay times (denoted by vertical bars in Fig. |3ft>)) also 
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FIG. 3: Measurements on a fluidized bed reveal two transi- 
tions: a transition in bulk properties at 4> g — 0.582 ± 0.004 
and an arrest transition at <p a = 0.594 ± 0.003 |2l|. (a) An 
increase near <j> g in the spatial correlation length £ of mobile 
regions (Fig. 2(a)), and a decrease in £ beyond <f> g to £ = at 
4> a - (£ is shown for AT=0.25 sec (•), 2.8 sec (o), and 50 sec 
(□).) (b) An increase in the average correlation time in bed 
images (Fig 2(b)). The inset shows a power law dependence 
of £ on t near 4> g . (c) An increase near <j> g in the decorre- 
lation time r DW s from the intensity autocorrelation function 
(1/e point) of transmitted laser light |20j|. The bold solid 
curve, the correlation time from (b) scaled by 3 x 10 s , be- 
gins to deviate from t D ws at (j> g . The inset shows that at the 
arrest transition, r DW s no longer increases, (d) A decrease 
beyond 4> g in the exponent a for the Mean Square Displace- 
ment (MSD) for intermediate times indicated by the dashed 
fine in the inset, which shows the time dependence of MSD. 
(e) An inflection in d<f)/dQ at 4> g and a divergence in dcj>/dQ 
as <f> a is approached, (f ) An increase in d(AP) jdQ beyond 4> g 
and a peak at <f> a , where d(AP)/dQ becomes independent of 
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increases as 4> increases; similar behavior has been found 
in a lattice model |24|. The slowing of the dynamics is 
similar to other glass forming systems Q. We find that 
the length and time scales for <j> < <p g are linked by a 
power law, as shown in the inset of Fig- Elb). A power 
law was also obtained in a study of a lattice model |24|, 
but with an exponent value 1/4 rather than the 3/2 given 
by our observations. In contrast to [24| the lengthscale 
plotted in Fig.Etb) is obtained for fixed AT = 0.25 sec. 

Additional signatures of glass formation are revealed 
by the bed behavior on short time scales. Measurements 
of the intensity-intensity correlation function g^ 2 '(t) of 
light multiply scattered as it travels through through the 
bed yield the time scale r DWS (Fig. Etc)), which in- 



creases until ■ 



with the same functional form as the 



decorrelation time of image pixel intensity r (Fig. ED 3 ))- 
The correspondence in the behavior for times differing 
by a factor of ~ 10 5 indicates that, as in other glass for- 
mers [lj, the macroscopic and microscopic dynamics be- 
have the same prior to the glass transition and become 
decoupled as the system forms a glass. 

Information on particle motion at small length scales 
was obtained using Diffusing Wave Spectroscopy (DWS) 
theory [2jj to invert the autocorrelation curves and ob- 
tain (Ar(t) 2 ), the mean square displacement (MSD) of 
the grains at times so short that particles moved only 
~ 0.01% of their diameter. For 4> <C <t>g, the short time 
dynamics are well described by (Ar(t) 2 ) oc t a with a w 2 
(Fig. Etd) inset), indicating that all grains undergo bal- 
listic motion between collisions |22t Erjj]. For (j> 0.54, 
as <j) approaches </> fl , the curves develop a region with 
a < 2 at short times, indicating that at these times the 
particles remain in contact with neighbors. However, at 
larger times, a w 2 (cf. inset of Fig. Etd)), indicating 
that the particles move ballistically. As <j> approaches 
(j) g , the MSD develops a plateau region at longer times, 
indicating the particles are caged; caging has been ob- 
served in hard sphere systems that approach the glass 
state 0, 0, 0|. However, the plateau observed in our 
MSD measurements indicates a caging length of order 
10 _4 <i, three orders of magnitude smaller than what is 
observed in these systems. For </> > </> ff , the slope at in- 
termediate times decreases, indicating that particles are 
no longer able to break from their cages, and a rapidly 
decreases (Fig. Etd)). For > cp a , all particles are im- 
mobile on macroscopic length scales (Fig Eta)- (c)). 

The bulk properties of the bed change at (f> g : d<fi/dQ 
reaches a minimum (Fig. Et e ))j an d the pressure drop 
AP across the bed begins to decrease rapidly, dropping 
well below the buoyant weight of the bed [27|. This is 
shown by the derivative d(AP)/dQ, which increases at 
4> g until the increase is arrested at tp a (Fig.Etf)), beyond 
which the bed has defluidized and AP becomes linear 
with Q, following Darcy's law [i^. We propose that the 
minimum in d<p/dQ and the increase of d(AP)/dQ be- 
yond 4> g are further signatures of glass formation: once 



the bed begins to behave like a glass close to 4> g , the ne- 
cessity for increasingly cooperative particle motion [Io| 
makes it unable to continue to pack sufficiently (increase 
cj)) in response to changes in Q as it would in the fluidized 
state where AP is maintained close to unity. Since the 
formation of a glass triggers the rapid drop in AP, which 
leads to the final macroscopic arrest of the bed at 4> a , we 
conclude that defluidization is a consequence of the hard 
sphere glass transition. This explains why <p a seen in 
our water fluidized bed and the gas fluidized bed exper- 
iments of Ojha et al. 01 is independent of particle size, 
container aspect ratio, and fluidizing medium. 
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FIG. 4: Microscopic motion in the glass state is halted by 
a slight increase in Q, jamming the system, (a) Q is de- 
creased (at 9.6 mL/min 2 ) to 12.5 mL/min, at which point it 
is increased by 1 mL/min (indicated by the rightmost dashed 
line). The glass state (when AP < 1) is achieved at Q — 24 
mL/min, indicated by the leftmost dashed line, (b) The 
speckle field, measured (30 ms exposure) along a line with 
each pixel corresponding to a single coherence area, becomes 
time independent after jamming occurs, (c) Hysteresis in <f> 
for increasing and decreasing Q after an increase in Q at Q « 8 
mL/min. The hysteresis disappears for cj> > cj> g . 



Jamming — Even in the arrested state, microscopic 
motion persists, as shown in the time evolution of images 
of a line of scattered light |2|| (Fig.^Jb)). The intensity 
of each pixel fluctuates even for <j> > <j> a , indicating that 
microscopic motion persists. However, for <f> > <j) g all 
microscopic motion can be jammed (stopped) by slightly 
increasing Q (for the conditions of our experiment, DWS 
can detect motion of only 1 nm for any scattering par- 
ticle). The jamming presumably establishes the stress 
backbone of the system [3(J and explains the hysteresis 
seen in Fig.QJc). Once the system is jammed, increases 
and decreases in Q below the onset of fluidization 0] do 
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not change 6. Such jamming: has been associated with 
glass states |3(J, |3l|, |32( and has been studied in fluidized 
beds H3. 

Conclusions — We have shown that a fluidized bed 
exhibits essential features of glass formation: rate de- 
pendence on final state, dynamical heterogeneity, rapid 
increase in time scale, and jamming. Further, the dy- 
namics of the hard sphere glass formation control the 
defluidization of the bed: beyond 4> g the system can- 
not pack sufficiently to accommodate changes in flow and 
must therefore arrest at <p a > 4> g with 4> a independent of 
particle size, container aspect ratio 18], and fluidizing 
medium. The volume fraction transition values depend 
weakly on surface properties |19j and strongly on cohesive 
effects ■ Finally, we speculate that the Random Loose 
Packed volume fraction, RLP ~ 0.56 [3^|, plays a role in 
the onset of heterogeneity, similar to the concept of onset 
transition observed prior to the glass transition [3j| . 

A fluidized bed is a simple system that allows fine con- 
trol, so it is an ideal system for studying glass and jam- 
ming transitions and for informing theory Further, 
understanding glass behavior in a fluidized bed can in- 
form fluidized bed des ign, which is important in many 
industrial applications |15j . 
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